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Abstract 



The well-known idea to construct domain wall type solutions of field equa- 
tions by means of an expansion in the width of the domain wall is reexam- 
ined. We observe that the problem involves singular perturbations. Hilbert- 
Chapman-Enskog method is used to construct a consistent perturbative ex- 
pansion. We obtain the solutions to the second order in the width without 
introducing an effective action for the domain wall. We find that zeros of the 
scalar field in general do not lie on a Nambu-Goto trajectory. As examples we 
consider cylindrical and spherical domain walls. We find that the spherical 
domain wall, in contradistinction to the cylindrical one, shows an effective 
rigidity. 
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1 Introduction 



Dynamics of topological defects is to be deduced from underlying nonlinear 
equations describing evolution of fields the defect is composed of. In general, 
it is a formidable task and one is forced to resort to approximate analytical 
solutions or to numerical solutions. The results are of great interest for 
particle physics (e.g. for dynamics of a flux-tube in QCD H), for field- 
theoretical cosmology (e.g. cosmic strings 0), and for condensed matter 
physics (e.g. domain walls in magnetics vortices in superconductors [|J 
or in superliquid Helium ||, defects in liquid crystals ||). 

In particular, in the context of relativistic field theoretical models there 
is a long-standing, interesting and difficult problem how to calculate evo- 
lution of vortices and of domain walls. Theoretical descriptions of vortices 
and domain walls show many similarities, and the two kinds of objects can 
be treated by essentially identical methods. The problem has become even 
more interesting after the seminal paper |7[] has appeared, in which it has 
been pointed out that the vortex dynamics can approximately be described 
in terms of a string model. Apart from purely numerical computations, see 
e.g. || H PH , the predominant approach is the effective action method, 
started in the paper [IT| and pursued in numerous papers till nowadays, see 
e.g. 0, |13|, 0, |T|, [16], [17], gj| [HJ. The idea is to reduce the complicated 



field-theoretical dynamics to much simpler dynamics of a string or a mem- 
brane by regarding some degrees of freedom as inessential. In the first step, 
these degrees of freedom are eliminated by the assumption that they have 
definite values determined from a conditional stationary action principle for 
the original field-theoretical action. The condition usually adopted is that 
zeros of the scalar field have apriori fixed, arbitrary trajectory in Minkowski 
space-time (identified with the world- volume of the membrane). Inserting 
these conditional solutions into the original action one obtains the effective 
action for the membrane. In the second step one has to find solutions of the 
equations of motion for the membrane, derived from the action principle for 
the effective action. Let us stress that in this approach calculation of the 
effective action is the necessary intermediate step. Without it, one has only 
the whole set of conditional solutions which exist for all trajectories of the 
membrane, including unphysical ones. It is just the role of the second step 



to pick physical trajectories for the membrane. See the paper [X3] for an 
explicit formulation and [pO] for a discussion of this method. 
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Another group of approaches aims at deriving directly from the original 
field equations dynamical equations for certain characteristics of the topo- 
logical defect, e.g. for an average radius in the case of a ring-like vortex pll], 



or position of zeros of the Higgs field [22]. In these approaches one does not 



have to calculate the effective action as the necessary intermediate step. To 
this group belongs also the approach based on a polynomial approximation 
to the fields constituting the defect |25| - in this approach also a width of 
the domain wall is a dynamical variable. 

In the present paper we would like to discuss one more method of cal- 
culating the evolution of domain walls or vortices directly from the original 
field equations: an expansion in the width of the domain wall or a vortex. 
This method has been tried in many papers, e.g. in |TT|-|T^|, always in con- 



junction with the effective action method. In the present paper we show that 
the expansion in the width itself, when appropriately carried out, is sufficient 
to determine evolution of the defect. No effective action is necessary. We 
will describe our formulation of this method in the case of the domain walls. 
Analogous calculation for vortices will be presented as a separate paper. 

Our main steps are moreless the same as in the previous papers: we 
introduce a special coordinate system co-moving with the domain wall and 
we formally expand the sought for solution of the field equation in the positive 
powers of the width of the domain wall. (The actual expansion parameters 
are dimensionless ratios of the width with curvature radia of the domain 
wall.) The first new feature specific for our formulation is that we do not 
tie the co-moving coordinate system with position of the zeros of the scalar 
field. It is tied to a Nambu-Goto membrane - this fact follows from the field 
equation as a consistency condition. However this membrane in general is 
not the locus of the zeros of the field. The position of the zeros is determined 
afterwards, from the explicitely constructed approximate solution of the field 
equation. 

Second, we observe that the expansion in the width is an example of so 



called singular perturbations |24J]. In such case one has to take into account 
certain consistency conditions for the perturbative expansion. Using Hilbert- 
Chapman-Enskog method for constructing perturbative expansion in the case 



of singular perturbations, see |25| for a lucid description of it, we obtain 
the domain wall solution directly by approximate solving the original field 
equation, without any need to introduce the effective action. 

The plan of our paper is as follows. In Section 2 we introduce the co- 
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moving coordinates and we present the field equation. This section is in- 
cluded mainly for convenience of the reader. In Section 3 we apply the 
Hilbert-Chapman-Enskog method to calculate the domain wall solution in 
the perturbative expansion in positive powers of the width. In Section 4 
we present simple examples: solutions describing cylindrical and spherical 
domain walls. Section 5 contains ending remarks. 



2 Transformation to the co-moving coordi- 
nates 

We shall consider domain walls in a well-known model defined by the follow- 
ing Lagrangian 

L = -- V ^<S>d^-^ 2 -—) 2 , (1) 

where (rj^) = diag(—l, 1, 1, 1), and A, M are positive constants. The corre- 
sponding field equation for the scalar field <3> has the form 

M 2 

^$-2A($ 2 -— )$ = 0. (2) 
4A 

There are two vacuum values of the field <£> equal to ±$ , where $ = 
M/2y/\. The domain wall is given by the field $ smoothly interpolating be- 
tween the two vacua. At each instant of time the field $ vanishes somewhere 
inside the domain wall. We assume that the locus of these zeros is a smooth 
surface S. We shall call it the core of the domain wall. Well-known example 
of the domain wall, with the core given by the (a; 1 , x 2 ) plane, is given by the 
following exact, static solution of Eq.(2) 

<$> = %tanh(^-), (3) 

where l = M^ 1 . The width of this domain wall is of the order l , and energy 
density is exponentially localised around the (x 1 ,^ 2 ) plane. 

For a generic domain wall, we parametrise the world- volume £ of the core 
(a 3-dimensional manifold embedded in Minkowski space-time, whose time 
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slices coincide with S) as follows 



(X») (r,a\a 2 ) 



( r \ 

X\r,a\a 2 ) 
X 2 (r,a\a 2 ) 

U 3 (r^V 2 ) J 



(4) 



where r coincides with the laboratory frame time a; , and a 1 , a 2 parametrise 
the core S at each instant of time. 

We shall also need another surface S attached to the domain wall. In 
general S ^ S. In the next section we shall obtain a Nambu-Goto equation 
for S, so this surface can be regarded as the Nambu-Goto type relativistic 
membrane. The surface S is related to the core S by a simple formula (51) 
derived in Section 4. In the following we shall call S the membrane co- 
moving with the domain wall. Let us stress that the membrane is merely an 
auxilliary mathematical constructs - the physical object is the domain wall. 

The world- volume of the membrane will be denoted by £, and we shall 
parametrise it as follows 



(r,a\a 2 ) 



I r \ 

Y\r,a\a 2 ) 
Y 2 (r,a\a 2 ) 
\V\r,a\a 2 ) ) 



(5) 



where r coincides with the laboratory frame time x°. The parameters a 1 , a 2 
in formula (5) are in principle independent of the ones in formula (4). We use 
the same notation for convenience, and also in anticipation of a one-to-one 
correspondence between S and S given by formula (51). 

As usual, we introduce a special coordinate system (r, a 1 , a 2 , £) co-moving 
with the domain wall. These coordinates are defined by the following formula 



» = Y^T,a\a 2 )+in»(T,a\a 2 ), 



(6) 



where x M are Cartesian, laboratory frame coordinates in Minkowski space- 
time, and (n fl ) is a normalised space-like four- vector, orthogonal to the E in 
the covariant sense, i.e. 



n 



,(T,a 1 ,a 2 )Y»(T,a 1 ,a 2 ) = 0, 



1, 



(7) 



4 



where a = corresponds to r; a — l,a — 2 correspond to ex 1 , a 2 ; Y^ = 
dY^/dr, etc. The four-vectors Y tT ,Y i<T i,Y >(T 2 are tangent to S. For points 
lying on the membrane £ = 0, and the parameter r coincides with the lab- 
frame time x°. For £ 7^ 0, r in general is not equal to the lab-frame time 
x°. Notice that definition (6) implies that £ is a Lorentz scalar. In the 
co-moving coordinates the membrane is described by the simple, Lorentz 
invariant condition £ = 0. 

The next step is to write Eq.(2) in the new coordinates. It is convenient 
to introduce extrinsic curvature coefficients K ab and induced metrics g ab on 
£: 

K ab = n,Y^ g ab = Y^ b , (8) 

where a, b = 0, 1, 2. The covariant metric tensor in the new coordinates has 
the following form 

(9) 

where a, (5 — 0, 1, 2, 3; a = 3 corresponds to the £ coordinate; and 

G ab = N ac g cd N db , N ac = g ac -^K ac . (10) 

Thus, G# = 1, Gi. a = (a=0,l,2, as in (7)). It follows from formula (10) 
that 

V=G = y/=gh(T,aW,Z), (11) 
where as usual g = det[g ab ], G = det[G a p\, and 

h(r,a\a 2 ,0 = 1 ~^K + ^(K^ - K b a K a b ) - ^ h abc e^ K a d K b e K c f . (12) 

Also g can depend on r, a 1 , a 2 . For raising and lowering the latin indices of 
the extrinsic curvature coefficients we use the induced metric tensors g ab , g ab . 
The inverse metric tensor G a/3 is given by 

(13) 
where 

Q ab = ( N -lyc gcd ( N -iyb_ (M) 



[G a p] - 



G ab 
1 



[G 



afl] 



Qab 
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Explicit formula for (N 1 ) ac has the following form: 

(N-'r = I {sT[l - ZK b b + \e{K h b Kj - K d b K h d )\ (15) 

+ i{l-iK b b )K ac + eK a d K dc } 

(this is just the matrix inverse to [iV a &]; by definition it has upper indices). 

In general, the coordinates (r, a 1 , a 2 , £) are defined locally, in a vicinity of 
the world- volume E. Roughly speaking, the allowed range of the £ coordinate 
is restricted by the smallest of the two curvature radia of the membrane in 
the local rest frame of a piece of the membrane. We assume that the domain 
wall is sufficiently narrow so that outside of the region of validity of the 
co-moving coordinates there are only exponential tails of the domain wall, 
i.e. the solution is exponentially close to one of the two vacuum solutions. 
Detailed discussion of the region of validity of the co-moving coordinates has 
been given in [ 23 ]. 

In the co-moving coordinates the field equation (2) has the following form 

-7=5^^°^ + > d ^ - 2A <* 2 - £>• = °< < 16) 

where u° = r, u 1 = a 1 , u 2 = a 2 . 

Let us rescale the scalar field and the £ coordinate 

$(^) = ^>K,s), t=^a, (17) 

where <j)(u a , s) and s are dimensionless. Then, equation (16) has the form 

2 1 » :( V=5G-|*) + i?$ + ^-^-l)# = 0. (18) 



M 2 y/=Gdu aK du b> 2ds 2 2hdsds 

We see that the expansion parameter, i.e. the width l = 1/M, appears in 
front of the highest order derivative with respect to r present in this equa- 
tion in such a manner that in the zeroth order approximation the r derivative 
drops out (as well as in the first order approximation in this particular ex- 
ample). It is well-known that in such cases straightforward perturbative 
expansion in positive powers of the expansion parameter in general is not 
correct p4|j . 
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The energy-momentum tensor in our model has the following components 

T» u = d»<$>d u <S> + TfL, (19) 

with L given by formula (1). The laboratory frame energy E is given by the 
integral 

E = J d 3 xT 00 . (20) 

In this formula the field $, which is a scalar with respect to coordinate 
transformations, can be regarded as a function of the (u a , £) variables. With 
the help of a formula for differentiation of composite functions T 00 can be 
written in the following form 

T oo = [(iv-y^ + n ^] 2 (21) 

+ 1 -g ab (N- 1 r(N~ 1 ) bd <!>^ 4 + i($ 5 ) 2 + V(*), 

where a = 0,1,2; n° is the /i = component of the four-vector (n M ); the 
potential 

\/($) = ^($ 2 -$ 2 ) 2 , (22) 

and the matrix A r_1 is given by formula (15). 

Also the volume element d 3 x is expressed by the co-moving coordinates, 

d 3 x = y/=jj h ^pQ M a l da 2 . (23) 

Observe that on the hyperplane of constant x° the r variable in general 
becomes a function of £; this is because 

r + £n°(r, a 1 ) = x° = constant. (24) 

In general, constant x° does not correspond to constant r, except for the 
membrane where £ = 0, hence x° = r. This complicates very much calcu- 
lation of the integral over £ because in general we do not know the explicit 
form of the r dependence - for this one would have to know explicit solutions 
of the Nambu-Goto equation for the membrane. 

There is a particular case in which the calculation of the energy is rela- 
tively simple: when each piece of the membrane associated with the domain 
wall is at rest at the initial time x°. Mathematically, this means that the 
first derivative of Y with respect to x° vanishes. Then, one can show that 
n° = 0, hence x° = r for all £. 
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3 Solution by the Hilbert-Chapman-Enskog 
method 



Following the earlier papers we expand the field (j> in the non-negative powers 
of 1/M, 

*^ (0) + > , + ^ <2) + ^ <3, + -. < 25 > 

and we solve equation (18) perturbatively. The expansion parameter is 1/M 
and not 1/M 2 because 1/M in the first power is present in the h and G ab 
functions after passing to the s variable, formulae (17). 

Inserting the perturbative Ansatz into equation (18) and expanding the 
l.h.s. of it in powers of 1/M we obtain the following sequence of equations 
for 0^: 

^ S V°) + 0(°)[1-(^)) 2 ]=O, (26) 
\d 2 s ^ + [1 - 3(0(°)) 2 ]0( 1 ) = Kd s ^\ (27) 

1^(2) + 3 (0(O))2] (2) = _ 2D (3) (O) (28) 

+2K b K£ sd s <f>W + 30W(0«) 2 + K a a d s ^\ 



1^0(3) + [! _ 3(0(o))2 ]0 (3) = AK a c K b a K c b s 2 d s ^ - 2n®0( 1 ) (29) 
+2K b a K a b sd s <f>W + 6^ U (1 U (2) + (0 {1) ) 3 + Kd s ^ 2 \ 
In these equations 

□< 3 > = -^L^-(^A) (3°) 

y^gdu a du b 

is the three-dimensional d'Alembertian on the world-volume £ of the mem- 
brane. Notice that only Eq.(26) is nonlinear. The other equations are linear 
inhomogeneous ordinary differential equations for (f)( n \n > 1. 

In general, the range of the variable s in these equations is restricted, 
4^£_ < s < 4f£+> because the coordinates (w a ,£) do not cover the whole 
space-time. The values of £± are found from the condition that h(u a ,C,) > 
for all £ in the interval (£+,£_). In general £± depend on u a . It is clear that 
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the function h with generic and g a b can vanish for finite £, see formula 
(12). From a mathematical point of view there is no problem - just one 
has to add more coordinate maps, smoothly matched with the coordinates 
(u a , £), in order to cover the whole space-time. From a practical, calculational 
point of view this would be a significant complication with rather little gain 
in understanding the physics of the not-too-much-curved domain wall we 
consider in this paper, because the additional coordinate maps would cover 
the uninteresting from physical point of view region in which the $ field is 
exponentially close to one of the vacuum values. For this reason we do not 
introduce the additional coordinate maps, and therefore we do not determine 
the field in the region where the domain wall exponentially merges into the 
vacuum. 

Then another problem arises, namely how to impose the boundary con- 
ditions for our solutions, because for finite £ = £ ± the field <E> in general 
does not reach the vacuum values. This problem can be circumvented in 
the following way. We shall find solutions of equations (26-29) assuming 
that s varies from — oo to +00. Notice that these equations formally make 
sense for all s. The solutions have physical meaning only in the restricted 
range of s. Because we presume that the domain wall field for £ close to 
£± is exponentially close to the vacuum values, we require that for large s 
(f)^ is exponentially close to ±1, and (fi n \n > 1, are exponentially close to 
zero. In this manner the boundary conditions at £ = £± are fixed uniquely, 
eventhough in the somewhat implicit manner. 

The zeroth order equation (26) has the well-known solutions 

<pf(s) = ±tanhs. (31) 

They obey the boundary conditions specified above. We take as the 
starting point for the approximation procedure. In the case of planar domain 
wall this coincides with the exact solution (3). The corrections given by 
cf)( n \n > 1, can be regarded as being due to the curvature of the domain 
wall. 

Let us pass to the first order equation (27). On the r.h.s. of it we have 
the function 

= ( 32 ) 

coshes 

The crucial step in the Chapman-Hilbert procedure is to observe that the 
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operator 

L = \dl + 1 - 3(0(°)) 2 , 

present on the l.h.s. of Eqs. (27-29) has a zero eigenvalue. The coresponding 
eigenfunction coincides with the function (32). (It is identical with a zero 
mode for a kink in the 1+1 dimensional version of the considered model (1).) 
Multiplying the both sides of Eq.(27) by the function (32) and integrating 
over s in the interval s G (— oo, +oo) we find that the integral on the l.h.s. 
vanishes, while the r.h.s. is equal to Therefore, we have to accept the 

consistency condition for our expansion 

K = 0. (33) 

It can be shown that this condition is equivalent to Nambu-Goto equation 

n (3) F M = o (34) 

for a relativistic membrane. This is the reason for calling the surface S 
the membrane. Solutions of the equation (33) for Y l have been discussed 



in numerous papers, see,e.g., |LJ], |2H|| . From Eq.(33) with fixed initial data 
one can calculate the Y % functions, and subsequently the extrinsic curvature 
coefficients K a b and the metric g a b- Then, from the condition h > one can 
determine £±. 

Because of condition (33) the r.h.s. of Eq.(27) in fact vanishes. Then, 
the general vanishing at the infinity solution of Eq.(27) has the form 



(35) 



cosh 2 s ' 

where C does not depend on s. 

The next equation to be solved is Eq.(28). Because (jr+ given by formula 
(31) does not depend on u a , and K% = 0, there are simplifications on the 
r.h.s. of Eq.(28): 

+ [i - 3 (,<»>) v > = + <;:;;;;. m 

Again on the l.h.s. we have the operator L, however this time the integration 
with the zero mode (32) gives the identity 0=0 because the r.h.s. of Eq.(36) 
is an odd function of s. 
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General solution of Eq.(36) can be found by a standard technique consist- 
ing of first finding two independent solutions of the homogeneous equation 
L0( 2 ) = 0, and next constructing a Green function, see e.g. \2B\. As the two 
linearly independent solutions of the homogeneous counterpart of Eq.(36) we 
take 

# = -4r, ( 37 ) 

coshes 

= -sinh{2s) + -tanks + - — S —. (38) 
8 8 8 cosh z s 

As the Green function we take 

G(s,x) = 2^\s)(j)?\x)G(s -x)- 2<$ ) {s)<t>f' \x)[Q{s - x) - Q(-x)}. 
The general solution of Eq.(36) has the form 

0( 2 ) = affi + + [ +CO dxG{s, x)f(x), (39) 

J — oo 

where a(u a ), f3{u a ) do not depend on s, and 

f = 2K b a KZ^— + 3C>^ (40) 
cosn z s cosn°s 

is the r.h.s. of Eq.(36). Formula (39) gives 

<pM = Cl ( w », s )0j 2 )( s ) + C2 (^, s )^\ s ) (41) 

where 

Cl (u a , s)=a- QC 2 h{s) - AK b a K a b h{s), (42) 

c 2 (u a , s)=(3- - 4K b a KZf 3 (s), (43) 

cosrrs 

The functions f%(s) are defined by the following formulae 

Ms) = f S dx^^$\x), f 2 (s) ee [ S dx^—^ix), (44) 
Jo cosh°x Jo cosn z x 

A « s -/>;dk- (45) 
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(2) 

The functions /i, /2 are odd , while /3 is even. The function 02 is given by 
formula (38). Because this function exponentially grows for s — > ±oo, the 
coefficient function c 2 has to vanish in this limit. Therefore, we have to put 

P = 0, (46) 

while a is still arbitrary. 

One more restriction comes from the third order equation (29). Integrat- 
ing it with the zero mode we obtain the following condition 

□< 3 >C + (^ - V)K h a Kt C + A c 3 = (^ - 1) KK\Kt. (47) 

Notice that the function a is not present in (47) - it drops out because 
an integral giving the coefficient in front of it vanishes. The coefficients in 
Eq.(47) has been obtained from computations of integrals over s. 

Condition (47) restricts possible choices of the function C{u a ). It has the 
form of a non-linear 2+1 dimensional wave equation. The function C can 
be regarded as a scalar field on the Nambu-Goto manifold determined from 
equation (33). 

Let us summarize results of our calculations. Up to the second order in 
1/M 

1 C la 
^ = tanhS+ M^s + JP^hKs (48) 

M l coshes cosh b s cosn z s 
where C obeys Eq.(47), while for a no explicit restrictions have been ob- 
tained till now. However, such restrictions might appear as the consistency 
condition following from the fourth order equation. Notice that Eq.(47) for 
the coefficient C introduced in the first order solution (35) has been obtained 
from the equation (29) two order higher. Thus, in order to claim that formula 
(48) gives the domain wall solution up to the order 1/M 2 one should also 
find 0( 3 ** and check the consistency condition coming from the fourth order 
equation. 

The fourth order equation has the following form 
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+ 34 O) (0«) 2 + 60i O) (1) (3) + 3(0«) 2 (2) (49) 

+2K b a K£ sd s (j)W + 4(K b a K«) 2 s 3 d s (j)^ + AK a c K b K c b 8%^. 

The r.h.s. of this equation has already been simplified by taking into account 
the condition (33) and the fact that 0+^ does not depend on u a . 

0(3) 

is given 

by a formula analogous to (39). Repeating once again the calculation leading 
to the consistency conditions we obtain from Eq.(49) the following condition 

U^ a + l -K b a K a b a-f^a = Q. (50) 

Thus, the fourth order equation (49) gives the restriction (50) for a. 

The final step in our procedure is to make a convenient choice of the 
initial position and velocities of points of the membrane. Its evolution at 
later times is then determined by the Nambu-Goto equation (33). Let us 
recall that the membrane is the auxilliary object merely, used to define the 
co-moving coordinate system, so that the choice of the initial data has no 
physical meaning. The most natural choice is that at the initial value of r 
the position of the membrane as well as its velocities coincide with those of 
the core. By the definition of the membrane, for its points s = 0. On the 
other hand, we find from formula (48) that for the core (i.e. for the points 
where <p vanishes) 

S - S "^-Ti-W + 0{M ^ < 51 » 

Therefore, our choice of the initial data for the membrane means that sq = 
0, d T so = at the initial value of r = r . This implies that 

dC 

C(u a ,r ) = 0, _( M Vo) = 0, (52) 
da 

a(u a ,r) = 0, — K,t o )=0. (53) 

Notice that the more general solution a = —MC is not acceptable, because 
one should not mix different orders of the perturbative expansion. With the 
initial conditions (53) equation (50) implies that 

a = (54) 
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for all r. On the other hand, the term on the r.h.s. of equation (47) acts as 
an external, curvature dependent force which in general leads to nontrivial 
r-dependence of the function C even for the trivial initial data (52). 

Thus, the domain wall solution constructed to the second order in 1/M 
with the Hilbert-Chapman-Enskog method has the form (48) with a = and 
C determined from the equation (47) with the initial data (52). 

The formula (51) gives the one-to-one mapping between the core and the 
membrane: the point {u a , s = 0) of the membrane is mapped to the point 
(u a ,s = so) of the core. 



4 Examples 



In this Section we specify our general solution to describe simple, symmetrical 
domain walls. 

(a) The planar domain wall. Let us start from the simplest domain wall: 
the planar, static one. In our formalism this means that the membrane is a 
motionless plane. In this case K ab = 0, so that the external force in Eq.(47) 
vanishes. Therefore C — 0, and the solution (48) coincides with the exact 
solution (3). Energy per unit square of the membrane is equal to 

i = -M*§, 
\S\ 3 °' 

where \S\ denotes the area of a piece of the membrane. 

(b) The cylindrical domain wall. Now the membrane is a cylinder of the 
radius r. We assume that the cylinder does not move as a whole. Using the 
formulae of the Section 2 one can easily calculate the extrinsic curvatures 
K ab and the metrics g a b- The Nambu-Goto equation has the following form 

f 1 

0, (55) 



where the dots denote derivatives with respect to r. This equation has the 
"energy" first integral given by the formula 



E 



N-G 
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where E G is a constant. We also find that 

2 2(E N ~ G ) 2 
K K = 5J = 1 r4 ' . ^°^X C = 0. (56) 

Thus, the external force in the equation (47) vanishes, and the membrane 
coincides with the core for all times - the zeros of the scalar field follow the 
Nambu-Goto trajectory. 

The solution cf) is given by formula (48) with C = a = and K^K h a as in 
the first of formulae (56). 

The energy for the cylindrical domain wall solution is also easily calcu- 
lated with the help of formulae (20), (21), (23). Because of the difficulty 
mentioned at the end of Section 2, we calculate the energy at the initial 
r = r , assuming that the core is then at instant rest. Thus, 



r(r ) = r(r ) = 0, 
and C has the initial value (52). At later time 

r(r) = r cos — , 
r 

as it follows from Eq.(55). It is easy to see that E N ~ G = r . 

Let us notice here that our formulae are applicable only for r/r G 
[0, 7r/2 — 5], with some 5 > 0. The reason is that for r/r approaching 
7r/2 it no longer true that on the inner side of the domain wall the field 
$ is exponentially close to the vacuum. The value of 5 can be estimated 



analogously as in the first of papers [23]. 



The result for the energy per unit square of the membrane is the same as 
for the planar domain wall. Thus, the extrinsic curvature of the cylinder does 
not influence the total energy per unit square of the membrane (to the order 
1/M 2 ). However, the density of energy across the domain wall is changed in 
comparison with the planar case. Its maximum still lies on the core, but it 
is flattened: the maximal value of the energy density is equal to 



2 ol M 2 rg j ' 
where /3(C)) = |/n2 — | ~ 0.30, and r is the initial radius. 
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(c)The spherical domain wall The membrane is a sphere of the radius r. 
The Nambu-Goto equation has the following form 



+ - = 0. (57) 



1 -r 1 

Its solutions have been considered in, e.g., WB, EH). r(r) is given by an elliptic 



function and it is oscillating with r. The "energy" integral is given by the 
formula 

e n-g _ r 



r 



;,2 



where E N ~ G is a constant. Again we consider evolution of the domain wall 
only in a part of the first quarter of the cycle of the oscillations predicted by 
Eq.(57), before the spherical domain wall collapses to the center. 
For K?K h „ we have the formula 



a 



K«K* 6 6 (^) 2 (r R) 



The external force in Eq.(47) does not vanish: 



KKKt = -—-^ = (59) 



r 



I _ ^2^3/2 



Formulae (58), (59) have been obtained by taking u 1 = <p,u 2 = 6, where 
<f), 9 are the usual azimuthal and polar angles parametrising the spherical 
membrane. The vector n has been directed to the outside of the sphere. 
Then s > (s < 0) corresponds to the outside (inside) of the spherical 
membrane. 

Equation (47) has now the following form 



At the initial instant r = tq : C = 0, r = tq, r = 0, and the acceleration 
C = — (tt 2 — 6)/?"o < 0. Thus, for r > C will become more and more 
negative. It follows from formula (51) that for the core s = Sq = —C/M > 0. 
This means that the core lags behind the Nambu-Goto membrane. 
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For the total energy per unit square of the membrane for the spherical 
domain wall we obtain a positive correction due to the curvature of the sphere 



As usual, we have calculated the energy at the initial instant; then r = r ,r = 
0; C = 0, d T C = 0. We see that the total energy E differs from the zeroth 
order result by the constant term 87rdo3?o/M. The maximum of the energy 
density at the initial instant lies on the core, and the maximal value of the 
energy density is equal to 



Thus, the energy density is more flattened than in the case of cylindrical 
domain wall. 

The facts that the core shrinks more slowly than the Nambu-Goto mem- 
brane, and that the energy per unit square has the positive correction in- 
creasing with the curvature, mean that the spherical domain wall exhibits an 
effective rigidity. 

5 Remarks 

(a) In our opinion, the main advantage of the method proposed in this paper 
is that it yields the domain wall solution without invoking the effective action. 
In this manner we avoid many problems of the effective action method, e.g. 
handling a model with higher derivatives. Another advantage is that the co- 
moving coordinate system is tied always to the Nambu-Goto type membrane. 
Equations of motion for this type of membrane are probably the simplest ones 
to be dealt with. 

In comparison with the polynomial approximation, [ffijjl , in which one 
approximates the field inside the domain wall by a simple polynomial, the 
present approach has the advantage that it is based on the expansion with 
respect to the clearly defined parameter, i.e. the width. Therefore, it is easy 




(61) 



where 
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to point out the limit in which the zeroth order solution is more and more 
accurate: C/M and K£K b a /M 2 tending to zero. The polynomial solution 
never is exact. 

(b) The corrections 4>^ l \ (f)^ to the zeroth order solution do not change 
significantly the picture of the evolution of the domain wall obtained from 
the zeroth order solution. For example, it is easy to see that in the case 
of spherical domain wall the difference between the radia of the core and 
of the membrane is of the order Iq/tq - it is smaller than the width of the 
domain wall. Nevertheless, it is nice to know how to construct the more 
accurate solution. Moreover, there are questions which one cannot answer 
without detailed knowledge of the internal structure of the domain wall. Let 
us mention here the problem of excitations of the domain wall. They can 
be regarded as bound states of certain particles with the domain wall. In 
order to calculate spectrum of such excitations one has to know details of 
the shape of the scalar field inside of the domain wall. 

(c) In the present paper we have completely neglected the possibility 
that the field $ can have a component oscillating with the frequency of the 
order M. Because M is the mass of the scalar particle present in the non- 
topological sector of the model (1), one should expect that such component is 
in general present. It could be regarded as describing an emission or absorp- 
tion of the scalar particles by the domain wall. Such oscillating component 
can be calculated in an extended perturbative expansion, which could be 
constructed, e.g., by generalising considerations of Section III of the first of 
papers 
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The problem of excitations of the domain wall and the problem of radia- 
tion probably deserve a separate investigation. 
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